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Employing neXus, one of the most recent simulation programs for heavy ion collisions, we investigate in detail 
the Hanbury-Brown Twiss correlation function of charged pions and kaons for central reactions 158 AGeV Pb-|- 
Pb. For this study we supplement the standard simulation program by electromagnetic interactions. We find 
a wide distribution of freeze out times which corresponds to a broad distribution of source sizes. Strong space- 
momentum space correlations as well as the electromagnetic interaction between the correlated pairs and the 
environment, which is of the same size as that between the hadrons of the correlated pair, modify strongly the 
correlation function. It can still be well parameterized by C(k) = l + \e~'^^ ^ , where k is the relative momentum. 
R is, however, not a simple function of the source radius but the complicated space-time structure of the source 
as well as the space-momentum space correlations and the final state Coulomb interactions are encoded. The true 
source radius, as given by the rms radius of the emission points, and R differ by up to 30%. 
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I. INTRODUCTION 



Which particle and energy densities are reached during 
an ultra relativistic heavy ion reaction? This is one of the 
present key questions in this field. If the density exceeds 
a critical value, theory predicts that a plasma, consisting 
of quarks and gluons, is produced whereas below that 
density the matter stays in its hadronic phase, consisting 
of mesons and baryons. 

Theoretical predictions on the density and the energy 
density obtained in present (SPS) or future (RHIC and 
LHC) heavy ion reactions are rather vague because the 
understanding of the initial phase in which projectile and 
target become decelerated is still a theoretical challenge 
without a reliable answer yet. Therefore the present ef- 
forts concentrate on an experimental determination of 
these densities. 

It is all but easy to measure coordinate space depen- 
dent quantities in heavy ion reactions. Single particle 
spectra contain direct information on the momentum 
space only. Therefore more complicated observables have 
to be employed to study the density. Probably the most 
powerful but certainly the most employed approach is 
based on a suggestion of Goldhaber, Goldhaber, Lee and 
Pais (GGLP) 1 1 1 who used the Hanbury-Brown and Twiss 
(HBT) effect [2]. It is based on the fact that the ampli- 
tudes of two indistinguishable processes interfere. If an 
object emits two identical bosons or fermions with mo- 
menta ki and k2 which are registered by two detectors 



the dependence of the cross section ^jfc^^tifc, on the rel- 
ative distance between the detectors or on the relative 
momentum of the registered particles is an image of the 
emitting source. This method has been proven to be 
quite powerful for measuring the size of stars [|| . 

In heavy ion physics the situation is much more diffi- 
cult: 

a) There is nothing like a static source which emits parti- 
cles. The system has a complicated geometrical structure 
and expands in space and time. Therefore not only one 
parameter (like in astrophysics the rms radius of a star) 
but the functional dependence of several parameters on 
time is encoded in the data. 

b) There is no unique definition of when a particle is 
emitted because the emission is in reality a disconnec- 
tion from the expanding system. Usually one assumes 
that after its last collision the particle can be considered 
as emitted. The space-time point of the last collision is 
called freeze out point. 

c) There are long range Coulomb forces which change the 
particle momenta after emission. 

Two approaches are possible in this situation: 
A) One neglects the complications, limits the problems 
with the geometry by using only particles which are emit- 
ted at midrapidity and parameterizes the source with a 
couple of parameters. The parameters are then deter- 
mined from the measured one and two particle cross sec- 
tions M. For a review on the present status of this ap- 
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proach we refer to the recent paper of Wiedemann and 
Heinz [||. The problem is that it is not easy to reveal 
the physical significance of the numerical values of these 
parameters. 

B) One employs numerical simulation programs which 
have been validated by reproducing the measured single 
particle spectra. These programs, however, are semiclas- 
sical in their nature and hence particles and not wave 
functions are propagated. Therefore one has to extend 
these programs by a quantum calculation (classical par- 
ticles do not interfere) if one would like to employ them 
for the prediction of the HBT correlation function. Here 
the rationale is as follows: If the extended program re- 
produces the correlation data, one can assume that it 
describes properly not only the momentum but also the 
coordinate space distribution of the particles. Hence the 
distribution of the particles in these simulation programs 
can be used to calculate densities and energy densities. 
The first suggestion to use the simulation programs to 
study the HBT effect has been advanced in Later 
Csorgo et al. |^ and Pratt et al. [Q have performed the 
first calculations using Fritjof, an event simulation pro- 
gram. In these calculations it has been assumed that a) 
the Coulomb interaction can simply be taken into ac- 
count by multiplying the two particle Wigner density 
by the Gamow correction factor and b) one can iden- 
tify the Wigner density of the nucleon with the classical 
phase space density. The latter is obtained by averag- 
ing the phase space points of the particles at the end of 
the simulation over many events. This has become the 
standard approach for a comparison of simulations and 
experiments . Already in ref. |^ it has been found that 
in heavy ion collisions the Gamow correction factor over- 
predicts the true effects and rather a classical trajectory 
calculation has to be employed. Recently several propo- 
sitions have been made to improve the treatment of the 
Coulomb interaction ||lc|-[T^ . The difference between the 
classical one body phase space density and the Wi gne r 
density of a emitted particle has been pointed out in . 

Here we study these both assumptions. In particular, 
using one of the up to date simulation programs which 
reproduces the single particle spectra quite well (section 
2), we study how the source parameters, obtained from 
the asymptotically observed correlation function, are re- 
lated to the true source. This includes a study of the 
source (section 3), a study of the Wigner density of the 
emitted particles and its relation to the classical phase 
space density (section 4), a study of the influence of the 
Coulomb fields (chapter 5), a study of the influence of 
the resonance decay and a study of the kaon correlation 
function. We will show that the relation between the 
source and its image (coded in the correlation function) 
is all but simple and that the source density as deter- 
mined by the correlation function may differ by a factor 
of two from that determined directly in the simulation 
program. 

Whereas the study of the difference between the 
Wigner density and the classical phase space density and 



the consequences for the interpretation of the correlation 
function in terms of a source radius is rigorous, the study 
of the influence of the Coulomb interaction can be only 
considered as a first step. It shows the order of mag- 
nitude of the corrections one has to expect but present 
computer power does not allow to supplement the simu- 
lation programs by the solution of a Schrodinger equation 
to calculate the (anti)symmetrized wave function of the 
correlated pions in the field of all other charges. 

II. THE METHOD 
A. neXus 

To study the space time structure of ultrarelativistic 
heavy ion collisions one has to rely on simulation pro- 
grams like VENUS neXus @ or URQMD 
These event generators describe the reaction from the ini- 
tial separation of projectile and target until the final state 
which is registered in the detectors. They assume that 
during the heavy ion reaction strings are formed which 
disintegrate into hadrons. The hadrons are treated semi- 
classically and can scatter among themselves elastically 
or inelastically before they are detected. For our study we 
use neXus (version 2.0 beta), the successor of VENUS, 
one of the most developed event generators, which de- 
scribes the single particle spectra of the observed mesons 
and baryons at CERN energies quite well. We refrain 
from describing the details of this model and refer to ref. 
p6[ . The program is used in its version without droplet 
formation in order to avoid the unknown Coulomb forces 
between droplets and hadrons. In using one of these mod- 
els we follow here the common believe that up to the last 
strong interaction collision (freeze out) the hadrons are 
well described by a semiclassical approach. 

B. Implementation of the electromagnetic 
interaction in neXus 

In order to study more rigorously the influence of the 
electromagnetic interaction on the two-particle correla- 
tion function we have to implement it directly into the 
simulation program. Because the average distance ^(^2) 
(where ^2 is the time when the later emitted hadron has 
its last collision) between the correlated hadrons is rather 
large (as we will see in part HI) it is sufficient to imple- 
ment the interaction on a semiclassical level [^,0 . Cor- 
related hadrons or pairs we call those pairs of particles 
which have a relative momentum smaller than 100 MeV. 
As we will see later the correlation function differs from 
one only for those pairs and hence only they carry infor- 
mation on the size of the system. 

The electromagnetic interaction of a particle i with 
another charged particle j is given by 
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dui 



(1) 



Here j runs over all charged particles of the system. The 
particle i has the mass to^, the charge and the 4- 
velocity Ui = {1, v}/-\/l — v^. 
by 



The tensor F-';' is given 



(2) 



where X'^ is the relative 4 - distance between the parti- 
cles i and j. 

The formula d) assumes point like particles and the 
absence of acceleration of the particle j. In our case 
this is a valid assumption because at high energies the 
electromagnetic force changes the momenta only little. 
Applying this formalism at each time step in the neXus 
program (which updates the particle positions and mo- 
menta at common light cone times) we create Coulomb 
trajectories of all charged particles. 

For the analysis presented later we store the mutual 
Coulomb force between all charged particles at all time 
steps. Whether two mesons form a correlated pair we 
know only at the end of the simulation. Using these 
stored forces we can separate the change of the rela- 
tive momentum of the pair particles kg due to their 
mutual interaction (Akpai^) from that due to their in- 
teraction with the environment (i.e. all other charged 
particles) (Ake„„iron)- The final relative momentum of 
the correlated pair \^ f inal is 
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(3) 



As we will see in this classical trajectory calculation the 
average change of the relative momentum of the corre- 
lated pair particles due to their mutual Coulomb interac- 
tion is about as large as that due to the interaction with 
the environment. This suggests that also in a quantal 
calculation both are equally important. The importance 
of the Coulomb interaction between a particle and its en- 
vironment can also be seen experimentally by cornparing 
the measured 7r+ and 7r~ spectra as done in figjlj. We 
see that at relative momenta smaller than 100 MeV the 
spectra differ considerably and that this difference is well 
reproduced by the neXus calculation including the classi- 
cal Coulomb interaction. This relative momentum has to 
be compared with the relative momentum k for which the 
correlation function C(k) becomes one and hence ceases 
to carry information on the source size. As we will see 
later k is as well around 100 MeV. Hence the influence 
of the Coulomb interaction on the correlation function 
cannot be neglected. 
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FIG. 1. 7r~ spectrum divided by the 7r+ spectrum, both 
obtained by the WA98 collaboration [11], as compared to the 
NeXus result. 



III. SPACE TIME STRUCTURE OF THE 
SOURCE 

The first question we address is when and where the 
particles, which are finally seen in the detectors, are cre- 
ated. We limit ourselves to pions emitted at midrapidity 
Vcm — 1 < y < Ucm + 1- For our study we use the reaction 
Pb -I- Pb at 158 AGeV for b < 3.2 fm. Fig.|top presents 
the distribution of freeze out times of charged mesons. 
We display this time separately for particles created by 
string decay, by resonance decay or inelastic collisions, 
and by elastic collisions. We see no sudden freeze out, 
rather a distribution which is almost flat between 3 and 
10 fm/c. Without rescattering we would have seen one 
peak at 1 fm/c (string breaking) and another one around 
5 fm/c {p decay ). 

The di stribution of transverse distances Rt = 
\/ -V of pions at freeze out with respect to the cen- 
ter of the reaction is displayed in the middle part of the 
figure. We observe as well a rather flat distribution fol- 
lowed by a very long tail which ranges up to 18 fm! Hence 
according to the simulation program rescattering made 
it impossible to define a unique radius of the system 
at freeze out which could hopefully be measured. The 
transverse momentum distribution of the pions does not 
strongly depend on their origin. This can be concluded 
from the bottom part of fig.|^, where the momentum dis- 
tribution of the pions is displayed. 

More interesting than the average properties of all par- 
ticles are that of those pions, which finally have a small 
relative momentum with respect to another pion of the 
same charge state. We will see later that the correlation 
function differs from 1 only for those pairs of particles 
which have a relative momentum smaller than 100 MeV. 
These pairs we call correlated pairs and only they carry 
information on the size of the system. 
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FIG. 2. Distribution of freeze out times, freeze out 
radii and freeze out momenta for charged tt's with 
ycm — 1 < y < ycm + 1 from neXus+Coulomb. Black squares 
show particles from string fragmentation, open circles show 
pions from elastic collisions and open triangles those from 
resonance decay or from inelastic collisions. 

The distribution of freeze out times ^2 of those pairs 
can be directly translated into a distribution of source 

sizes Rsourceit2) = j^T^aii prons ■ This dis- 

tribution of source sizes is displayed in figj| for two cases: 
for all pairs and for those pairs where none of the pions 
come from resonance decay. We see a broad distribution 
with a maximum at 8 (lO)fm and a mean value of 11.3 
(12.4)fm. For calculating the mean value we have omit- 
ted all particles beyond Rsource = 30 fm, i.e. pions which 
come from long living resonances like w or A. 

Thus there is neither a common freeze out time nor a 
common freeze out volume. 




FIG. 3. Distribution of source sizes for pairs with a relative 
momentum of less than 100 MeV. 

In fig.^ we compare the distribution of the relative dis- 
tance r{t2) at freeze out for the correlated pairs in com- 
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FIG. 4. Distribution of relative distances at freeze out for 
pairs with a final momentum of fc < lOOMeV in comparison 
with that of all pairs of positively charged pions. 

parison with that for all pion pairs. First of all, the aver- 
age distance is above 10 fm. We observe furthermore that 
the average value for the small relative momentum pairs 
(■\/ < r^{t2) > = 10.8 fm) differs considerably from that 
for all pairs (\/ < r'^{t2) > = 15.6 fm). Consequently, the 
pairs with a small relative momentum are not representa- 
tive for all pairs. This one has to keep in mind if one likes 
to interpret the correlation function in terms of physical 
source radii. 
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FIG. 5. Momentum-position correlation in the neXus 
model. 

The difference between the distribution of relative dis- 
tances of pairs with k < lOOMeV and that of all pairs is 
a consequence of strong space - momentum space corre- 
lations, in beam (z) direction as well as in the direction 
of the impact parameter (x). These correlations of freeze 
out momenta and positions are displayed in fig. ^ and 
are well known. They give rise to the identity of rapidity 
and space time rapidity. We see that even in the small 
rapidity interval around midrapidity, in which we per- 
form our investigation, these correlations are by far not 



negligible. 

Now, after we have described the source, we can inves- 
tigate to which extend its properties can be recovered by 
measuring correlation functions. 



IV. THE CORRELATION FUNCTION WITHOUT 
ELECTROMAGNETIC INTERACTIONS 

A. The formalism 

All semiclassical simulation programs propagate clas- 
sical particles which have a sharp momentum kyiand a 
sharp position ta- This is conceptually necessary in or- 
der to determine the sequence of collisions and the center 
of mass energy of each collision. On the other hand, two 
particles emitted with sharp momenta and sharp posi- 
tions do not interfere because their trajectories are dis- 
tinguishable. In order to model physics which is based on 
the interference of wave functions, like the HBT effect, 
one has to treat the particles as quantum particles for 
which the uncertainty principle is fulfilled. Therefore a 
transition from classical to quantal physics is necessary. 
In view of the construction of the simulation programs 
this transition can only be made after the last collision. 

This transition from classical to quantal physics re- 
quires more than the knowledge of k^iand ta provided 
by the simulation programs. In order to construct a wave 
function one needs at least one additional parameter (the 
width of the wave function) which is not provided by the 
simulation program. A popular choice for the wave func- 
tion is a Gaussian superposition of plane waves. This 
wave function requires only one unknown parameter, the 
variance L/4, and there is even hope that its physical 
relevance can be understood, as we will see. 

Two identical particles emitted at ti and at the 
source points ta and tb with the momenta centered 
around \iA and k^ observed in detectors with the mo- 
menta ki and k2 can have two indistinguable trajectories. 
Either the particle emitted at arrives with momentum 
ki (and that from tb with momentum k2) 



fAip,) = (ii)3/4,-(Pl-k^) = V4. 
ZTT 



(7) 



= y dpidp2 < ki,i2|ii,Pi >< Pi,ti\^A > 

• < k2,<4|t3,P2 >< P2,i3|$B > (4) 

or that emitted at arrives with momentum k2 and 
correspondingly that emitted at with ki 

A2 = J dpidp2 < ki,i2|ii,Pi >< Pi,ti|$B > 

• < k2,i4|i3,P2 >< P2,t3|*A > • (5) 

For the wave function of the emitted particle 

< pi,h\^A /A(pi)e'Pi^-4e-«Pi=ti/2™ (6) 
we assume a Gaussian form 



If L is large the momentum of the particle is close to 
the source momentum but we loose the information on 
the source position. If L is small we have the opposite 
scenario. 

< ki,i2|ti,Pi >= '5(3)(ki - p^)e-*Pi'(*2-ti)/2m is the 
propagator in momentum space. Because both trajec- 
tories are indistinguishable their amplitudes have to be 
added to obtain the correlation function 
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JdK{Al + Al) 
g{k - ks) + g(k + ks) + 25(k)g(ks)cos(2kr5) 



(8) 



(9) 



gik-ks)+g{k + ks) 

where k= 1^1^ ,K = (ki+k2) ,k5 = ^^^^ , = 

(r^ - Tb) and g{k) = {^)^^*e~'^ ^. In our calculation 
we update the coordinates of the first emitted particle 
until the emission of the second particle of the corre- 
lated pair (at ^2). Therefore cos(fc^tr'') = coskr where 
= {0, r(i2)} and the time component of the 4-vector 
of the relative distance disappears from the correlation 
function. If the two bosons have the same source momen- 
tum [kA—kB) and are measured in the detector with the 
same momentum we obtain an enhancement of C^^-'(k) 
by a factor of two as compared to distinguishable parti- 
cles. However, if the source momenta are different, the 
result is more complicated. 

If L is sufhciently small, the difference in the argu- 
ments of the exponential functions can be neglected and 
we obtain 



.g(k - ks)(l + e-^'^cos{2krs)) 



This approximation is called smoothness assumption. It 
is frequently used in the literature for the following rea- 
son: The measured correlation function can be well de- 
scribed by a Gaussian function iQ 

C(k) = 1 + Aea;p(-4i?2 . k^) (n) 

where A allows for the correction of a possible coherence 
in the source or for pairs of non identical particles. As- 
suming that the emission points have as well a Gaussian 
distribution in coordinate and momentum space 



S{ks,rs)o,e-^s'''e-''s'/^, 



(12) 



where kg and rs are the relative distances in momentum 

f 21 

and coordinate space, the convolution of Cg^^(k) with 
the source distribution function S(k5,r5) yields 

CsM(k) = 

/ drsdksS{ks,rs)g{k - ks)(l + e-^''^cos{2krs) 



= 1 



JdvsdksS{ks,rs)g{k-ks) 

,-k'(L+C) 



(13) 
(14) 
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The apparent source size AR^ is, in the Gaussian source 
approximation, four times the sum of the variance of the 
source C/4 and of that of the wave function squared L/4. 
In this approximation we find the foUowing relation be- 
tween the root mean square radius of the freeze out points 
\/ < f (^2)^ > and the apparent source radius 



< r{t2f >= = 3/4(i + C). 



(15) 



Thus in the smoothness approximation one has only to 
fit the experimentally observed correlation function by a 
function of the form of eq.O and can then use eq. |l^ 
to determine the desired < r(i2)^ >. Consequently one 
can avoid to perform a Fourier transform, a difficult task 
in view of the large experimental error bars. In addition 
this equation is easy to interpret. 

In the simulation programs we do not have a continu- 
ous source distribution. Rather we have to sum over all 
emission points: 



C(k) 



j:jdK\A{ + Ai, 



E/dK(|AlP 
which yields in the smoothness approximation 



(16) 





FIG. 6. 2 different sets of single particle wave func- 
tions squared | $i(r) |^ which give the same one 
body Wigner density J F^{r,k)dk. On the right hand side 
\Mr) f={J F^„-{r,k)dk)/5. 



well. The argument goes as follows: From the classical 
n-body phase space density we can construct the classical 
one body phase space density (we use the nonrelativistic 
formulas here because there the argument becomes more 
transparent) 



CsA/(k) = 



cos{2krss,t)) * 5(k - kg^^,) 



(17) 



pairs is the number of pairs in the event s. 
Now it is important to realize which information is 
needed to calculate a correlation function and which in- 
formation is to our disposal in these simulation programs. 
As seen in eqs. ^,p p^Jl7| we need 

a) the single particle wave function for the two correlated 
particles emitted from emission points with relative co- 
ordinates rs and kg. 

b) the distribution of emission points. Because in the 
derivation of the correlation function the source is 
treated classically these emission points r5,ks are ei- 
ther given by a source function like eq.|l^ or as a sum 
of delta functions in coordinate and momentum space as 
in eq.|l|,|l]. 

The simulation programs treat the rescattering of the 
hadrons like the scattering of classical billiard balls. Only 
initially, when the hadrons are created by string decay, 
quantum mechanical constraints are taken into account. 
Consequently, the simulation programs present the time 
evolution of a classical chaotic system whose initial con- 
figuration is inspired by quantum mechanics. Hence we 
have to our disposal the classical n-body phase space den- 
sity. Assuming that emission takes place at freeze out the 
classical n-body phase space density allows to calculate 
the emission points. 

It has been argued in the past [ pTjp^ that the knowl- 
edge of the classical n-body phase space density is suf- 
ficient to determine the single particle wave function as 



F^,{vM = J26'-'\k-k,it)) 5(3)(r_r,(t)). (18) 



Averaged over many events |£T| or over one event |2^] 
the classical one body density can be identified with the 
quantal one body Wigner density 



(r , k, i) = ^ g(k, r, k. (t) , r, (t) ) 

PI) 



(19) 



where ^(k, r, ki{t),Yi{t)) is the Wigner density of the sin- 
gle particle wave function $(k, ki(t), ri(t)) defined as 



9{k,vMt)Mt))= J ^e'q^$*(k + q/2,k,(0,r,(i)) 

(20) 



(27r)3 

$(k-q/2,k,(t),r,(t)). 



Now one assumes that all ki{t),ri{t) are identical. Then 
we can identify 



F^{r,k,t)=n-g{k,r) 



(21) 



Knowing g(k, r) one can calculate the amplitudes Ai 
(eqg) and A2 (eq|§. 

There are, however, three shortcomings in this argu- 
ment: 

1) There is an infinity of single particle wave functions 
<I'(e q.po[) which give the same one body Wigner density 
(eqjljyr They yield, however, quite different correlation 
functions. Fig.^jshows as example two sets of single parti- 
cle wave functions which give the same one body Wigner 
density but different two body correlations C(k). Thus 
eq. O is one out of an infinity of choices of a single 
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particle wave function for a given one body Wigner den- 
sity. One may wonder if this is a good choice because it 
ignores all correlations between the ki(i)'s and ri(i)'s. In 
addition it is difficult to motivate why until freeze out the 
particles are treated as point like in coordinate and mo- 
mentum space whereas thereafter they are well described 
by a wave function for which AfcAr » h. It is more nat- 
ural to replace the precise momenta and positions of the 
particles by the "most classical" wave functions around 
the classical ki(t)'s and ri(t)'s, i.e. Gaussians which ful- 
fill the uncertainty principle. In general, the ignorance 
of L and the impossibility to derive from the known one 
body Wigner density the single particle wave functions 
are two ways to express the same unknown physics. 

Averaging over many events raises even more ques- 
tions: 

2) Two particles emitted next to each other (i.e. those 
which are finally interfering strongly) have a mutual 
Coulomb repulsion, which introduces a two body correla- 
tion between the pair even on the classical level. Averag- 
ing over many events (which means to allow that parti- 
cles from different events interfere) destroys this correla- 
tion. Thus event averaging destroys correlations already 
present on the classical level. 

3) There is no proof that the classical n-body phase space 
density averaged over many events is equivalent to the 
quantal one body Wigner density. On the contrary there 
are good reasons that this is not the case, if the size of the 
system is comparable to the width of the single particle 
wave function: If all particles are sitting on top of each 
other a classical n-body density is point like whereas the 
Wigner density has to obey the uncertainty principle. 

For the determination of the correlation function from 
experimental data the denominator is obtained by event 
mixing. For our theoretical studies this is not necessary. 
Therefore we would not like to comment here on the con- 
sequences of the event mixing procedure for the correla- 
tion function. Numerically we observe quite a difference 
in the results if we replace the denominator of eqs,16 1^ 
by the corresponding quantity for mixed events. This is 
caused by mutual Coulomb interactions (see 2)) which 
get lost by event mixing. 



B. The Correlation Function in neXus without 
electromagnetic interactions 

As discussed, for a Gaussian source (eq.[l^) with no 
space-momentum space correlations there is a simple re- 
lation between the true source radius and the apparent 
source radius (eq.p^. 

Is this observation also true for the source as given by 
the simulation programs? This is the seminal question 
which we will discuss in this and in the next paragraphs. 
In the simulation programs we are confronted with three 
different source radii 



• the c lassical source radius 

Rciass — \/ < ''(^2)^ > as given by the simulation 
program where < . . > means averaging over all pi- 
ous which are part of a pair. Their radius is taken 
at the time when the later emitted pair particle has 
its last collision. The distribution of r(i2)^ (i.e. of 
the freeze out points) is displayed in fig. ^ We 
see that this distribution deviates from a Gaussian 
form. It can be considered approximately as a sum 
of two Gaussians. 



• the true source radius Rtrue given by the rms 
radius of the system after the freeze out points 
have been convoluted with a Gaussian wave func- 
tion squared of a variance of L/4. This source ra- 
dius is called true because it is the relevant source 
radius if one wants to determine densities or en- 
ergy densities and it is also the source radius the 
correlation function should "measure". For L=0 

Rtrue Rciass- 

• the apparent source radius R as given by eqjTl|. 
This apparent source radius can be calculated from 
the correlation function using the formalism de- 
rived in the last section. This is the only radius 
which can be measured experimentally. 



collision -I- resonances, Mecn i0.51 " frr 
collision . Mean 1 1 .34^ f 




250 300 



FIG. 7. Distribution of the freeze out points. 

We start the investigation by calculating the correla- 
tion function (eq.^ or eqjl^) of identically charged pions 
at midrapidity for different values of the variance L/4. 
We discard for the beginning those pairs which contain 
pions coming from resonance decay. These correlation 
functions have a form which can be well described by 
eq. [l[ We can therefore determine R by fitting the sim- 
ulation results with a function of the form of eq.^. In 
figj^ we display this dependence of the fit parameter 
on L. We separate the results obtained with (eq.|l7|) and 
without (eq.l6) the smoothness assumption. 

Per definition for L=0 both agree. For larger values 
they differ by not more than 10% and hence the smooth- 
ness assumption is acceptable for the analysis of the re- 
sults of the simulation programs. For the unrealistic case 
that Vl is large as compared to the source we see the 
expected (see eq.|l5|) linear behavior, but this is not true 
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anymore if -s/X < Rdass, for the smoothness assumption 
as well as for the exact calculation. 



O Solution(eq 14) 

Approximate Solution (eq 15) 
Smootiiness assumption 



with x-p correlation 




O Solution(eq 14) 

Approximate Solution (eq 15) 
Smoothness assumption 




FIG. 8. The square of the apparent source radius R as 
a function of the variance L/4 of the single particle wave 
function squared. We display calculations with and without 
smoothness assumption. In the bottom figure we have arti- 
ficially decorrelated the positions and momenta of the parti- 
cles, is obtained by fitting the correlation function with 
C(k) = 1 + Xexp{-4k^R^). 

In order to understand the origin of this deviation from 
the expected behavior we decorrelate position and mo- 
mentum of the pions. We attribute to each pion with 
a given source momentum a position randomly chosen 
among the source positions of the other pions. The re- 
sult is displayed in the bottom panel of fig.^. We see 
now the expected behavior: The apparent source size in- 
creases if the width of the wavefunction and hence the 
true source radius increases. 

In fig. ^ we plot the true source radius as a function 
of the apparent source radius. As explained, both can 
be obtained independently in the simulation program. 
Without space-momentum space correlations the appar- 
ent source radius follows the true source size, the convo- 
lution of the classical source with the single particle wave 
function squared, which can be independently calculated 
in the simulation program. The absolute values are, how- 
ever, different. The reason for this difference is that the 
distribution of emission points (fig. |^) is not a Gaussian 
(despite of the fact that the correlation function can be 



well described by a Gaussian). 

If we include the correlations generated in the simula- 
tion program the apparent source radius is, however, a 
quite complicated function the true source radius. As- 
suming that the apparent source size is the true source 
size one overpredicts the density of the source by up to 
a factor of 2. Thus we conclude that for the reaction 
investigated here the true source radius cannot even ap- 
proximately be inferred from the apparent radius with- 
out a detailed knowledge of the space-momentum space 
correlations . We just like to mention that the experi- 
ments yield apparent source sizes which are considered 
as too small in order to be the true source size, a fact 
which, according to our results, may have its explana- 
tion in space-momentum space correlations. 
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• Smoothness assumption 
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FIG. 9. The true source size Rtme (as given by the sim- 
ulation program) as a function of the apparent source size 
R determined by fitting the calculated correlation function 
by C(k) = 1 + Xexp(-4k^R^) including (see fig. | top) and 
excluding (see fig. H bottom) space-momentum space corre- 
lations. 

We would like to add two remarks: 

a) The space-momentum space correlation in the sim- 
ulation program is not created by flow. It is already 
caused by the string dynamics and is modified later due 
to rescattering. Thus adding only a radial (or in plane) 
flow to models which have otherwise no correlations be- 
tween space and momentum space would not reduce the 
difference between R and Rtme- 

b) Being forced to introduce L as a free parameter we 
have lost the predictive power of the simulation programs 
as far as two body correlations are concerned. Instead of 
predicting the two body correlation function we can only 
state whether there is a value of L for which the simula- 
tion program reproduces the measured apparent source 
radius R. (As said the different values of Rtme are ob- 
tained by folding the same classical source distribution 
with the square of wave functions with different variances 
L/4.) 
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V. SIMULATIONS INCLUDING 
ELECTROMAGNETIC INTERACTIONS 

A. The momentum change due to the 
electromagnetic interaction 

After freeze out and during their way to the detectors 
the mesons change their momentum due to the Coulomb 
interaction. In fig. 10 we display the distribution of the 




70 80 
AP (MeV/c) 

FIG. 10. Distribution of the change of the relative momen- 
tum between the particles of a correlated pair after freeze out. 
We separate the change due to the electromagnetic interac- 
tion with the environment (black circle), with the other pion 
of the correlated pair (open cross), and the sum of both (open 
circle) . 

change of the relative momentum of the correlated pairs 
between t2 and detection, separated into the different 
origins. The average momentum change is of the order 
of 25 MeV and the interaction between the correlated 
pair contributes as much as the interaction of the pair 
particles with the environment. Even if these numbers 
are based on classical trajectory calculations, presented 
in section II, it is evident that calculations which take 
into account only one of these effects are not suited for 
the situation at hand. 



of the particle wave function 4'^{r) , where r is the rela- 
tive distance of the pair particles, onto the Coulomb wave 
function $^.^^'(r). 



C(k,K) = J d\\ < nTJ\r X r|0^ > 



$^y2'(r) is the solution of 



(22) 



(23) 



Assuming identical freeze out points for the two par- 
ticles \(j)^{r)\'^ = S^{r) we obtain the so called Gamow 
correction factor 



2ttt] 



exp[2TTri] — 1 



(24) 



with ij = jjLa/k , ji the reduced mass. Often the 
"Coulomb corrected" correlation function at freeze out 
is calculated by dividing the measured correlation func- 
tion by the Gamow correction factor. 

This approach has been recently critically studied and 
extended to a system of n bosons by Alt et al. . They 
have shown that in a certain limit the n-body wave func- 
tion can be approximated by the product of the relative 
wave functions of all pairs and that for small source sizes 
(« Ifm) this product can be replaced by that of the 
relative wave functions at zero relative distance. Then 
the Coulomb correction is nothing else than a product of 
Gamow correction factors. However, already for sources 
of 5 fm this product of correction factors overestimates 
the true Coulomb correction In view of the large 

source size we see in neXus (fig. 3) these results elucidate 
why the Gamow correction does not work well for heavy 
ion reactions at ultrarelativistic energies, as has been ob- 
served in 



B. Standard treatment of electromagnetic 
interactions 

As said already, the correlation function C(k) dif- 
fers from one only for pairs with a relative momentum 
k = \J (ki — k2)'^ < lOOAfeT^. Only these pairs carry 
information on the size of the system. The change of 
the relative momentum of a correlated pair due to the 
Coulomb interaction is of the same order of magnitude 
and therefore not negligible. 

This is well known since long and several attempts have 
been made to correct the two body correlation function 
for the electromagnetic interaction |]l0|-|l3||. The correla- 
tion function C(k, K) for a pair of particles emitted with 
a center of mass momentum K and an asymptotic rela- 
tive momentum k is given by the square of the projection 



C. The correlation function in presence of 
electromagnetic fields 

The above mentioned quantal approaches for the 
Coulomb correction are implicitly based on the assump- 
tion that the correlated pair particles are emitted very 
close to each other. Only then the environment changes 
the momentum of the pair particles by the same amount, 
leaving the relative momentum k unchanged. In this case 
one can solve eq.^ This is not the situation we see in the 
simulations. There the change of the relative momentum 
of a correlated pair particles due to the mutual Coulomb 
interaction is of the same order of magnitude as that due 
to the interaction with all the other particles. Then eq. 

is not anymore a good approximation even more it is 
not even calculable: the asymptotic relative momenta k, 
necessary to calculate the Coulomb wave functions, are 
only known at the very end of the simulation, after the 
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mutual Coulomb interactions have ceased, and are based 
in our approach on classical trajectory calculations. Thus 
we are left with the fact, that if quantal effects are impor- 
tant, our asymptotic relative momenta, based on a clas- 
sical calculation are not correct, on the contrary, if they 
are correct quantal effects are negligible and we can live 
with calculating Coulomb corrections classically, as done 
here. Thus the quantal Coulomb correction proposed in 
are incompatible with semiclassical simulation pro- 
grams. 

It has been proposed as well to replace in the Coulomb 
wave function the asymptotic momentum by the momen- 
tum at freeze out. For many applications this is certainly 
a very good approximation, but in our context it is not 
useful: this ansatz would eliminate the effect we would 
like to study: How the correlation function is modified by 
the interaction of the pair particles with the environment. 

We would like to mention another problem which ren- 
ders the application of a static quantal Coulomb correc- 
tion in the context of semiclassical simulation programs 
quite difficult: Coulomb interactions and strong interac- 
tions cannot be separated in time. The earlier emitted 
pair particle has moved already in a (strong) Coulomb 
field when the later emitted one has its last collision. 

Due to these problems even the most modest quan- 
tal approach, to describe the correlated pair by a two 
body wave function in which all other charges are treated 
in the Born Oppenheimer approximation, i.e. by sup- 
plementing in the Hamiltonian of eq. ^(r) by the 
Coulomb interaction with all other charged particles 
Sfc(^'^(''ifc' + ^'"(I'jfeiOj ^^^d then to solve the time 
dependent version of this equation, is not feasible with 
present day computers. 

What can we do in this situation? Because the ex- 
act calculation is not possible we employ approximative 
methods to study the influence of the electromagnetic 
interaction on the correlation function. These approxi- 
mations are crude, too crude probably to yield quantita- 
tively reliable results, but sufficiently precise to demon- 
strate that more precise methods have to be developed 
before a quantitative relation between the apparent ra- 
dius and the true source radius of the system can be 
established. 

Our approximation consist in assuming that the mo- 
mentum transfer due to the electromagnetic interaction 
is instantaneous. Under this assumption we study two 
cases. First we assume that the change of the momentum 
due to the electromagnetic interaction happens immedi- 
ately at freeze out. This result is then compared with 
calculations which assume that the instantaneous mo- 
mentum change happens at t, the time where in ncXus 
on the average half of the momentum change has hap- 
pened, t is of the order of 20 fm/c. Because the mo- 
mentum change due to the electromagnetic interaction 
is moderate it does not completely destroy the correla- 
tion function obtained by neglecting the electromagnetic 
interaction (as for example collisions with large momen- 
tum transfer would do). Rather one expects a smooth 



modification. 

The first case is easy to study. We have only to re- 
place ks by ks+q where q= Akpa„ -I- Akenvironment is 
the momentum change due to the electromagnetic inter- 
action. The second case requires a more careful study. 
As said, in order to create a correlation function one has 
to describe the mesons after freeze out by a wave func- 
tion with a finite width in coordinate and momentum 
space. Therefore, if we take the quantum character of 
the mesons seriously, we have now 4 different amplitudes 
which are visualized in fig. 9. 

Si = y" dpidp2dp3dp4dp5dpe 

< ki,i2|ii,P5 > <5(p3 + qi -ps) < P3,ii|0,pi > 

< Pl,0|$A >< k2,t2|il,P6 > <^(P4 + q2 - Pe) 

< P4,tl|0,P2 >< P2,0|$B > 

= < ki - qi,ti|$A >< k2 -q2, tile's > 

B2 = J dpidp2dp3dp4dp5dpe 

< ki,t2\ti,pe > (5(p4 + q2 - Pe) < P4,ti|0,pi > 

< Pl,0|$A >< k2,t2|<l,P5 > ^(P3 +qi -P5) 

< P3,tl|0,P2 >< P2,0|$B > 

= < ki - q2,ti\<^A >< k2 -qi,<i|$s > 
S3 = y dpidp2dp3dp4dp5dpe 

< ki,t2|<l,P6 > ^(P4 + q2 - Pe) < P4,tl|0,P2 > 

< P2,0|$B >< k2,t2|il,P5 > '5(P3 + qi - Ps) 

< P3,tl|0,pi >< Pi,0|$A > 

= < ki - q2, tile's >< k2 - qi,ii|$A > 

B4 = y dpidp2dp3dp4dp5dpe 

< ki,t2\ti,p5 > (5(p3 + qi -P5) < P3,ii|0,P2 > 

< P2,0|$B >< k2,t2\ti,pe > (5(p4 + q2 - Pe) 

< P4,tl|0,pi >< Pi,0|$A > 

= < ki - qi,ti|$s >< k2 - q2,ti\^A > ■ 

Calculating | S/p we find after integration over the 
center of mass motion: 

Re{BiB^) = Dcos(2q(r5 + 2kt/m))e-(^-^sf^-ce^ 
Re{BiBl) = Dcos(2k(r5 + 2qt/m))e-(q+ks)'i-k'L 
Re[BiBX) = Dcos{2Ys[k - q))e-(k-q)'i-k5'i 
Re{B2Bl) = Dcos{2vs{k + ci)))e-(^+'^)^^-^s''L 
Re{B2Bl) = Dcos(2k(r5 - 2qt/TO))e-(q-'*s)'i-k'L 
ReiBsBl) ^ Dcos(2q(rs - 2kt/m))e-*+^s)'i-q'i 
BiBr = 
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BaSI = De-*+q-ks)=L 
B3B* = i5e-*-q+ks)^L 
BiBl = i5e-(k+q+ks)=L 




kl k2 

FIG. 11. The 4 indistinguishable different amplitudes if the 
electromagnetic interaction happens instantaneously, with a 
momentum transfer of q2 and qi . 



where D contains the normalization. The exact result 
^(k) = ^|g'|2 is too complex to discuss the physics. 

Therefore we apply the smoothness assumption, i.e we 
assume that all exponential functions can be approxi- 
mated by e^*^'*^*l^^'5'^^. Then we obtain 

s=l i=l 

(1 + 0.5 X e-'-^-^^^^ cos(2r5,^,(k - q)) + 

0.5 X e-*+q)'L cos(2r5,,,(k + q)) + 
1 2 

cos(2rss,ik)cos(4qkt/m) + 
e~*^ ^ cos(2r5g_jq)cos(4qkt/m)) 

/(E E e-(k-q-ks..)^.). (25) 

If k is equal the correlation function is not anymore 
two but oc l+e^*l cos (2rq). For finite k and if q is small 

we find CsaO^) = 2(1 + e^^ ^cos(2kr5)) and hence up 
to a normalization constant the same result as in section 
IV. It is important to realize that the change of the mo- 
mentum due to the electromagnetic forces modifies the 
result in a nonlinear way. 

This general result will be compared later with that 
obtained under the assumption that qi is the momentum 
transfer of the particle emitted with the momentum k^i 
(and correspondingly that the particle emitted with k^ 
suffers a momentum transfer of q2). 

D. The correlation function in neXus including 
electromagnetic interactions 

Before discussing our results in detail we compare the 
calculated correlation function for 7r"'"7r"'" and 7r~7r~ cor- 
relations in fig.^. Here we have assumed that the change 
of the momentum due to the electromagnetic interaction 
happens at ^2 and we choose L — Ofm^. As observed 
experimentally (fig. |l|) the correlation function is almost 
identical for tt"*" and tt" pairs even including the Coulomb 
interaction in distinction to the single particle spectra 
(fig. |l]). Please note that (as discussed in section IV 
A) the correlation function would be strongly modified 
if one uses mixed events for the normalization. In true 
simulation events the correlated pair particles suffer from 
a mutual Coulomb interaction before freeze out, whereas 
in mixed events this not the case. 
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FIG. 12. Correlation functions for tt+tt"*" and 7r~7r~ pairs 
emitted at midrapidity in 158 GeV/N Pb+Pb central colli- 
sions 

We present now our results for three different scenar- 
ios. All are based on the assumption that the momentum 
transfer due to the Coulomb potential is instantaneous. 
In the first scenario we assume that the change of the mo- 
mentum due to the electromagnetic interaction happens 
at t2, in the second and third scenario that it happens 
at the time when in the classical trajectory calculation 
on the average half of the momentum change due to the 
electromagnetic interaction has occurred. The latter sce- 
nario is calculated for two different assumptions. First 
we assume that the particle emitted at A suffers a mo- 
mentum change of qi (and that emitted at B changed its 
momentum by q2). Then only the trajectories Bi and -B3 
of fig.|ll| are allowed. Finally we admit all 4 amplitudes 
Bi — S4. 

Fig.|l3| displays the correlation function under the as- 
sumption that the momentum transfer happens immedi- 
ately at freeze out. We separate 4 different source mo- 
menta to discriminate the different origins of momentum 
change: a)k5 momentum at freeze out (No Coulomb) 
b)ks -f Akenviron (Environment), the change of the rel- 
ative momentum of the pair due to the interaction with 
the environment. c)k5 + Akpair (Pairs), the change of 
the relative momentum of the pair particles due to their 
mutual Coulomb interaction and d)kfi„ai=ks + Akpair 
+ Akenviron (AH). The Correlation functions are rather 
similar and the fit yields very similar results for R. Due to 
the fact that we don't mix events for the determination 
of the denominator of eq.[l6|, this result was expected. 
Event mixing would create a shift of the maximum to a 
finite value of Qmt due to the lack of Coulomb repulsion 
between the particles from different events. 

If we assume that the momentum transfer happens at 
10 fm/c resp. 20 fm/c after t2 and admitting only the 
amplitudes Bi and B^ we observe the dependence of the 
true source on the apparent source radius R as displayed 
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FIG. 13. Correlation functions for 4 different final mo- 
menta (see text). We assume that the momentum transfer 
due to electromagnetic interaction happens at freeze out. L 
is assumed to be = 10/m^. 

m fig.|l|. For a given source size Rtme the Coulomb 
interaction increases the measurable apparent source size. 
The later the momentum transfer happens the larger is 
this change. For a given Rtme a momentum transfer at 
20 fm/c may increase of the apparent source size radius 
by as much as 25%. No Coulomb is identical with the 
curve displayed in fig.^ We like to mention that the 
change of the apparent source radius due to the Coulomb 
interaction is of the order of magnitude expected from the 
schematic study in ref. ^ . 

If we allow for all 4 different amplitudes Bi — B4 we see 
another time a quite different behavior. The correlation 
function has not anymore a Gaussian form as can be seen 
in figjl^. Fitting only the low momentum component 
by a Gaussian function one obtains R values which are 
slightly lower than if only the amplitudes Bi and B3 are 
allowed. The error bars of the fit are, however, large. 
We would like to mention that the non Gaussian form 
of the correlation function becomes only evident in our 
approach. If the correlation function is normalized by 
event mixing an overall normalization appears. Because 
the correlation function is normalized to 1 at large values 
of k this overall normalization absorbs the large k offset 
seen in fig.^Sj. 

The influence of the pions from resonance decay on 
the correlation function is rather weak although the res- 
onances can decay very late, when the distance to its 
pair partner is very large. This large distance increases 
Rtrue tremendously. Therefore we have excluded them 
up to now from the analysis. The apparent source ra- 
dius of all pairs as compared to pairs of particle where 
no pion comes from resonance decay is displayed in fig.[l6| 
as a function of L (L/4 being the variance of the wave 
function). This result is understandable: at the late time 
when the resonances disintegrate there is no other par- 
ticle close by and hence the value of qr(i2) is large. On 
the average this contribution is therefore negligible. 
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FIG. 15. Correlation function obtained for the general 
case where all amplitudes B\ — B4 (see text) are admitted. 
We assume that the momentum transfer takes place at t= 20 
fm/c. The correlation function is not anymore Gaussian. 
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FIG. 14. The true source radius Rtme (as given by the 
simulation program) as a function of the apparent source ra- 
dius R determined by C(k) = 1 + Xexp{—4k^ R^) including 
the Coulomb interaction (see text). We admit here only the 
amplitudes Bi and Bs (see text) and assume that the mo- 
mentum transfer due to the electromagnetic interaction takes 
place at t= 10 fm/c (left) or at 20 fm/c (right). 
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FIG. 16. The apparent source radius R determined by 
C(k) = 1-1- Xexp{—4k^ R^) as a function of L including and 
excluding the pairs in which at least one pion comes from 
resonance decay 
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VI. CORRELATION FUNCTION OF THE KAONS 

So far we have investigated pions. In the simulation 
programs we have sufficient to allow for the calcula- 
tion of a correlation function. The result is displayed in 
figjl^. We observe a very similar behavior as for the pions 
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FIG. 17. The true source radius Rtrue (as given by the 
simulation program) for kaons as a function of the apparent 
source radius R determined by C(k) = 1 + Aea;p(— 4k^i?^). 

because the space-momentum space correlations are sim- 
ilar. The kaon source is smaller than that of the pions: 
Rtrue as well as R are smaller. The cross sections 
are lower and therefore the freeze out happens earlier. 
Consequently the smaller source radius is a natural con- 
sequence of known physics. We observe as well a large 
discrepancy between Rtrue and R. 

VII. CONCLUSION 

We studied in detail two problems which one encoun- 
ters if one tries to interpret HBT correlation functions 
obtained from numerical simulations of heavy ion reac- 
tions at CERN energies: The influence of space - mo- 
mentum space correlations and the influence of final state 
Coulomb interactions. For this purpose we supplemented 
one of the presently available programs, neXus, by elec- 
tromagnetic interactions. 

We find, first of all, that there is nothing like a static 
source which emits particles. The particles come from 
sources of quite different sizes because the system ex- 
pands as a whole. Momenta and positions of the parti- 
cles are strongly correlated and hence the source is not 
chaotic. Both observations render the task rather diffi- 
cult to extract from the correlation function useful in- 
formation on the size of the system or on densities. The 
most one can expect to obtain are a time averaged values. 
However, even these time averaged quantities are neither 
unique (different particle species give different values) nor 
representative (pairs of particles with small relative mo- 
menta (to which the correlation function is sensitive) are 
emitted from a smaller source then arbitrary pairs). 



In order to calculate a correlation function in these 
semiclassical simulation programs one has to introduce 
(at least) one free parameter. Here we assume that the 
particles have a Gaussian wave function after their last 
collisions. Consequently, the simulation programs can- 
not predict the two body correlation functions. They 
can only determine for which variance the calculated two 
body correlation functions agrees with experiment. 

Knowing the variance of the wave function and the 
freeze out points we can calculate in the simulation pro- 
grams directly the average rms radius of the source. It 
is called Rtrue- AH correlation functions we obtain have 
the form C(k) = 1-1- Ae^"*^ , the same form as seen 
in the experiments. Thus the crucial problem is to relate 
the apparent source radius R with physically meaningful 
quantities like Rtrue- 

The simulation program tells us that there is no sim- 
ple relation between R and Rtrue due to four different 
reasons: 

1) The distribution of the emission points is not Gaus- 
sian, whereas (due to the the finite error bars in the simu- 
lation as well as in the experiments) the correlation func- 
tion is well described by a Gaussian and finer details of 
the source cannot be resolved in the correlation function. 
Therefore there is a difference between R and Rtrue al- 
ready of the order of 20%. Furthermore, the distribution 
of emission points and the distribution of sources (as dis- 
played in fig. 3) differ as well considerably. The latter is 
the relevant quantity if one would like to extract densities 
or energy densities. 

2) For the reaction investigated the strong space- 
momentum space correlations modify the relation be- 
tween R and Rtrue in an important way. Thus R does not 
measured the true size of the source but underestimates 
it by « 16% and hence the true average density of the 
system is 60% lower than that extract from R. This may 
explain the fact that the experimental source sizes, when 
analyzed under the assumption that the source is (al- 
most) chaotic are incomprehensibly small. The relation 
between R and Rtrue depends in a complicated way on 
the dynamics of the system and is not due to quantities 
like flow. Thus it is not evident how the information on 
space-momentum space correlations can be incorporated 
in more phenomenological approaches and, consequently, 
how this type of models can extract the true source radius 
from the observed apparent source radius. 

3) The Coulomb interaction changes the relation be- 
tween R and Rtrue as well. Assuming that our calcu- 
lation gives quantitatively reliable results, the difference 
between R and Rtrue can reach 30% and hence the den- 
sities differ by a factor of 2. For the influence of the 
Coulomb interaction we presented here only very approx- 
imative results. They are based on classical Coulomb tra- 
jectory calculations. They show that the Coulomb inter- 
action between the correlated pairs and the other charges 
is as important as the Coulomb interaction between the 
particles of the correlated pair and that the standard pro- 
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cedure to use Gamow correction factors is not justified, 
because the average particle distance at freeze out is too 
large. 

4) The difference between R and Rtme depends on the 
particle species because their cross section are different. 

Thus this study shows that two problems have to be 
solved before the correlation function measurements can 
provide useful information on the density of the system. 
We have a) to understand how the apparent source ra- 
dius is related to the time averaged true source radius and 
b) how the true source radius can be related to physical 
quantities like energy densities or densities. Here we have 
addressed the first question only and found that the dif- 
ference between R and Rtme is all but negligible. Hence 
a naive interpretation of the observed apparent source ra- 
dius as the radius of the system at freeze out is certainly 
premature. 

The Coulomb interaction produces still another effect 
if one employs the standard method of event mixing be- 
cause the absence of the Coulomb interaction in mixed 
events modifies the correlation function. This effect we 
have observed in the simulations but not discussed here. 
It will further add to the difficulty to relate the apparent 
source radius to the true source size. 
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